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Abstract. In SUSY models which violate R-parity, there exist trilinear lepton number violating (LNV) operators which 
can lead to neutrino masses. If these operators are defined at the unification scale, the renormalization group flow becomes 
important and generally leads to one neutrino mass much heavier than the others. We study, in a minimal supergravity 
(mSUGRA) set-up with two trilinear LNV operators and three charged lepton mixing angles, numerically how these 
parameters may be arranged to be compatible with neutrino oscillation data, and discuss some phenomenological observations. 
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MSUGRA MODELS WITH LEPTON 
NUMBER VIOLATION 

In the minimal supersymmetric extension to the Stan- 
dard Model (SM), a general superpotential contains si- 
multaneously lepton and baryon number violating terms 
[1]. Their presence leads to fast proton decay [2], so at 
least one of the above set of terms need to be forbidden, 
for instance by some discrete symmetry [3]. The most 
widely studied symmetry is R-parity [4], for which both 
sets of terms are projected out. An alternative, known as 
baryon triality [5], also suppresses fast proton decay but 
allows the set of lepton number violating (LNV) oper- 
ators. The presence of these LNV terms in turn leads to 
interesting phenomenological consequences such as neu- 
trino masses and neutrinoless double beta decay. It also 
opens up a large region of parameter space not avail- 
able in R-parity conserving models, because the light- 
est supersymmetric particle (LSP) in the latter needs to 
be electrically and color neutral to avoid cosmological 
problems [6]. For these reasons, it is important to study 
models with LNV. 

The renormalisable LNV superpotential is given by 

^LNV - ^XijkLiLjEk + X!j^LiQjDk^^iLiH,.il) 

We shall focus on the dimensionless, trilinear opera- 
tors Xjjii and Xlii^. The supersymmetry (SUSY) break- 
ing sector is specified using the minimal supergravity 
(mSUGRA) [7] assumption. We wiU attempt to account 
for the neutrino oscillation data [8] (see Table 1) by two 
LNV parameters. They are defined in a weak interaction 
basis in which the charged lepton Yukawa matrix Ye is 
symmetric but not diagonal, and the rotation to the di- 
agonal basis can be characterized by three angles. The 
parameter set that defines our model is therefore 

mo,Mi/2,Ao,sgn;U, @ Mx, 



TABLE 1. Neutrino oscillation data obtained in [8]. 

Am2^ = 7.9+[]-^^ ■ lO^^eV^ |Am^i | = 2.6±0.2 ■ lO^^eV^ 

sin20i2 = 0.31 ±0.02 sin^dij = 0.47+j]J]^ 

^■,„2a _n+0.008 
sin Wi3 — 0_oQ 



tanjS @ Mz, (2) 

where Ai,A2 G {A,j<., A/^^.}, andMx is the scale at which 
the electroweak gauge couplings unify, mo, /2 and Aq 
are the universal scalar mass, gaugino mass and SUSY 
breaking trilinear scalar coupling at Mx respectively, and 
sgn/i is the sign of jj,. The charged lepton mixing angles 
0''s are defined in the standard parameterization [9], 
and tanjS is the ratio of the higgs vacuum expectation 
values Vu/vci- It is important to note that because the 
PMNS [10] mixing angles are large, the charged lepton 
mixing angles are expected to be of (?{\). This means 
that despite the assumption of having two dominant LNV 
parameters at the weak interaction basis, after rotating to 
a diagonal charged lepton basis there will be many LNV 
parameters of similar order of magnitude. 



Neutrino masses 

Since lepton number is violated, the 3 neutrinos mix 
with the 4 neutralinos. The 7x7 neutrino-neutralino 
mass matrix ./^^ can be written as 

= "M, (3) 

where and m contain the lepton number violating 
contributions, and is the neutralino mass matrix. 
An effective 3x3 neutrino mass matrix can be obtained 




FIGURE 2. A mass insertion diagram which represents the 
radiative correction to my by ^HjRt and Xj^K. . 



by an 'electroweak' see-saw mechanism, in which the 
neutrahnos act as the see-saw. It is given by 

■^eff ^ my—m^^Qm^. (4) 

Renormalization and radiative effects 

Note that at tree level, my =03x3- The effective mass 
matrix is proportional to [11] 

{^^ff)ij - A,-A;, 

Aj = fivi~v,ifii, i^{e,fi,T}, (5) 

so depends only on the sneutrino vevs v,- and /i; 

bilinear parameters. Even though they are not present 
at Mx, they will be generated upon renormalization to 
lower scales. A pseudo Feynman diagram for the gener- 
ation of fii is shown in Fig. 1. 

The rank one structure of ^^^^ means there is only 
one massive neutrino at tree level. The inclusion of ra- 
diative corrections is necessary to provide a second mass 
scale required by the measured neutrino mass squares 
differences. The one loop corrections to {.^n}:] [12] is 
given by j{d^N + where 

i5^N)i,j = {T.D)ij-{.^N)ik{'^L)kj, (6) 

and So and are mass corrections and wavefunction 
renormalization respectively. A mass insertion diagram 
for the radiative correction of niy is shown in Fig. 2. 

While the bilinear parameters which control the tree 
level mass scale are generated via renormalization group 
running, the trilinear parameters already present at Mx 
contribute to the radiative correction of .^n. Due to the 
large logarithmic factor coming from renormalization 
from Mx to ~ 1 TeV, the tree level contribution turns 



out to be dominant. The ratio m'"" /niy"'^ can be approx- 
imated to be [13] 

mT _ _^ aGC/rln^(Mx/Mz) 
m'v""" ~ 107!;Mi/2(Ao - AttanjS ) 

where ric is 1 for Xijk and 3 for A/^^, / is a dimensionless 
function of ^(10)[14], and is mass matrix of scalar 

/. A rough scan in the mSUGRA parameter space gives 
ratios of i??'(30 — 200). On the other hand, consistency 
with neutrino oscillation data requires the mass ratio of 
the two heaviest neutrinos to be at most of ^(5). This 
implies for our simple model to work, the combined 
contribution to the dynamical generation of the bilinear 
parameters must be suppressed to decrease the mass ratio 
described above. The charged lepton mixing angles also 
need to be fixed. The best fit parameters are obtained by 
a numerical procedure discussed next. 

NUMERICAL PROCEDURE 

We restrict ourselves to the SPSla [15] point 

Mi/2 = 250GeV, mo = lOOGeV, Ao = -100GeV, 
sgnjU = +1, tanjS = 10, 

(8) 

and vary only the 2 LNV parameters and the 3 charged 
lepton mixing angles 0'. After specifying these param- 
eters at Mx, the Lagrangian is rotated to a diagonal 
charged lepton basis using the 0''s. A modified version 
of SOFTSUSY [16] code, including the full set of 1- 
loop LNV MSSM RGEs contiibutions, is used to run to 

^(500)GeV at which the EWSB conditions ai-e im- 
posed, and the neutrino mass matrix is calculated. All 1- 
loop corrections to the latter are calculated numerically 
except for the CP-even and CP-odd neutral scalar contri- 
butions, which are computed in an analytic expansion in 
the relevant LNV parameters to avoid fluctuations from 
strong numerical cancellations. The optimization of the 
best fit parameters is performed using MINUIT with the 
neutrino oscillation data displayed in Table 1 . 

In Table 2, a selected set of best fit parameters is dis- 
played. We note that no attempt to exhaust all combi- 
nations of X and X' couplings is made. There are other 
parameter sets which fits the neutrino oscillation data, 
but are not presented here since they violate other exper- 
imental constraints [17], such as leptonic FCNCs. 

In the table, we display values of and fine tuning 
measure Aft of the best fit points. The tuning is much 
more severe for an inverted hierarchy because in the for- 
mer case, instead of having a symmetry to relate the two 



TABLE 2. Set of best fit parameters at SPSla, assuming 

both the normal and inverted hierarchies. See text for more 
discussions. 
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1^33 = -2.50- 10" 
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0.460 


0.389 


0.305 


8.09 
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4.07 


10- 


-5 


1.990 


1.082 


0.632 


8.10 




^533 = -3.41 • 10" 
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9.87 


10- 


-5 


0.448 


0.400 


2.891 


12.4 




A;22 = -1.1410" 
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4.06 


10- 


-5 


1.193 


0.191 


1.174 


11.8 




A;22 = -8.98 10- 


5 






1.03 


10- 


-4 


2.107 


0.175 


1.181 


9.44 




A,;22 = -8.60 10- 
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Mii 




4.10 


10- 


-5 


0.998 


0.282 


0.418 


8.00 




Inverted hierarchy 
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A2 
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eis 


Aft 




A533 = -5.69 10- 
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■^233 




1.36 


10- 


-4 


1.558 


0.815 


0.146 


755 


0.01 


X533 = -5.69 10- 
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Mil 




1.32 


10- 


-4 


1.388 


0.760 


0.141 


758 


0.06 


Aa3 = -5.68 10" 
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Al23 
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-4 
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-0.758 
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726 


0.05 
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Mil 




1.25 


10- 


-4 


0.135 


0.102 


0.798 


988 


0.43 


A;22 = -1.94 10" 
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1.47 


10- 


-4 


3.032 


0.087 


0.932 


743 


2.85 


A;22 = -1.96 10" 
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M23 




1.46 


10- 


-4 


0.144 


0.094 


0.690 


736 


0.52 



neutiino masses to be (almost) equal, the LNV parame- 
ters need to be arranged to make the tree level mass scale 
quasi-degenerate with the scale of radiative corrections 
in ^gff The lack of an underlying flavour symmetry 
also means the observed near tri-bi maximal mixing [18] 
should be regarded as accidental. 

A few general observations on possible LHC signa- 
tures are in order. The LNV parameters are typically too 
small to affect decay of the sparticles produced, except 
for the LSP. On the other hand, the charged lep- 
ton mixing angles means that the LSP is expected to de- 
cay in many different channels with similar branching 
ratios. For the SPSla spectrum we examine, the LSP can 
decay with a displaced vertex of ^(0.1 )mm. Such dis- 
placements should not be immediately obvious, but may 
be searched for to provide support for this type of mod- 
els. However if the LSP is a stau for instance, it should 
undergo 2 body decay at the interaction point. 

SUMMARY 

We have discussed SUSY with lepton number viola- 
tion as an alternative to neutrino masses. In high scale 
constructions, renormalization effect typically leads to 
neutrino masses incompatible with neutrino oscillation 
data, unless the trilinear LNV parameters are arranged to 
suppress the dynamical generation of the bilinear LNV 
parameters responsible for the tree level my. We per- 
formed a numerical study to fit the oscillation data in an 
mSUGRA model, assuming two dominant triUnear LNV 
couplings defined in a weak interaction basis, and discuss 
coUider effects that could be observed at the LHC. 
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